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Abstract 

We analyze the capability of charm production in deep-inelastic ep scattering 
at HERA to constrain the gluon distribution g(y, /i 2 ) of the proton. The de- 
pendence of the theoretical predictions for the charm structure function F£ on 
the mass factorization scale fi and the charm mass is investigated, i 7 ^ seems 
to be well suited for a rather clean and local gluon measurement at small 
momentum fractions y. 
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We analyze the capability of charm production in deep-inelastic ep scattering at 
HERA to constrain the gluon distribution g(y, [i 2 ) of the proton. The dependence 
of the theoretical predictions for the charm structure function f„ c on the mass 
factorization scale fi and the charm mass is investigated. F£ seems to be well suited 
for a rather clean and local gluon measurement at small momentum fractions y. 



Recent measurements of the structure function _F 2 ep at HERA rather directly 
determine the (light) quark densities in the proton at very small momentum 
fractions y. On the other hand, extractions of the gluon distribution from F% 
scaling violations are more indirect and require assumptions on the evolution 
kernels at small y, e.g. of the validity of the usual NLO evolution equations, 
which themselves are a matter of investigation presently. An obvious candidate 
for a more direct gluon determination is the longitudinal structure function Fl , 
but in practice this quantity is not easily measured sufficiently accurate over 
a wide kinematical range. In this context it is interesting to recall that F2 
itself is expected to contain an appreciable part directly sensitive to the gluon 
density at small Bjorken-x, namely its charm production contribution F^a. In 
the following we discuss some theoretical aspects of this observable relevant to 
the determination of g{y,jjL 2 ); for the experimental status at HERA see ref. □. 

In next-to-leading order (NLO) perturbative QCD, the electromagH|etic 
(exchange of one photon with virtuality Q 2 ) structure function F£ readscl 
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Here m c (e c ) denotes the charm mass (charge), and as is the strong coupling 
constant. The renormalization scale has been put equal to the (MS) mass 
factorization scale fi. The lower limit of integration over the fractional initial- 
parton momentum y is given by y mm = ax = (1 + Am 2 /Q 2 )x, corresponding to 
the threshold s = Am 2 of the partonic center-of-mass (cm.) energy squared. 
The NLO coefficient functions and c^ 1 ' have been calculated, in ref. 0, and 
a convenient parametrization of these results has been provideda in terms of 
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In leading order (LO), 0(as), F 2 C is sensitive only to g(y,/i 2 ) via the well- 
known Bcthc-Hcitler process j*g — > ecu. A comparison of the various con- 
tributions to F.f in NLO shows that for the physically reasonable scales /z, 
(X ~ 2m c . . . \/Q 2 + 4mj? (see below), the quark contribution in (1) — which is 
not necessarily positive due to mass factorization — is very small, about 5% 
or less. Therefore F£ does represent a clean gluonic observable also in NLO. 
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Figure 1: The ^-dependence of the coefficient functions ej ^ 1 ', derived in refs. 3—5 for two 
values of Q 2 with m c = 1.5 GeV. The scheme-dependent quantity c 2 g is given in the MS scheme. 

The ^-dependence of the gluonic coefficient functions 3-5 is recalled in 
Fig. 1 for two values of Q 2 typical for deep-inelastic small- a; studies at HERA, 
as in the following using m c = 1.5 GeV. The comparison of the NLO coefficients 
c- 1 ' and c' 1 ' with the LO result reveals that potentially large corrections 
arise from regions where is small, i.e. from very small and large partonic 
cm. energies. These corrections are due to initial-state-gluon bremsstrahlung 
and the Coulomb singularity at small s, and due to the flavour excitation (FE) 
process at s > Am 2 (rj 3> lA. For a more detailed discussion, including the 
quark contributions, see ref. 13. Large FE-logarithms have been resummed — 
at the expense of losing the full small- s information of (1) — by introducing 
a charm parton density, leading to the so-called variable-flavour schemeB. For 
another approach to high s see ref. 0. The importance of these corrections in 
the HERA small- a; regime considered here will be investigated below. 

The first question to be addressed in order to judge the phenomenological 
usefulness of F 2 C as a gluon constraint is whether or not the available NLO 
expression (1) is sufficient for obtaining results which are stable under variation 
of the (unphysical) mass factorization scale. It has been argued U that one 
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should use jj, ~ 2m c , since /i is supposed to be controlled by s and the integrand 
in (1) is maximal close to the lower limit, s ~ 4m 2 . The range of significant 



contributions in s, however, broadens consider, 
below), hence /i ~ \7Q 2 + 4m^, chosen in refs 
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with increasing Q 2 (see 
appears at least equally 



reasonable. Therefore we estimate the theoretical uncertainty of F 2 C in NLO 
by varying the scale between /j = (2)m c and /i = 2\/Q 2 + 4m;?, see Fig. 2. 
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Figure 2: The dependence of F£ on the scale fi for m c < fi< 2y/Q 2 +4m% at selected values 
of x and Q 2 . GRv| and CTEQ20 (at x = 10 -3 , 0.04) parton densities have been employed. 

One finds that at small-x, x<!0~ 2 , the scale variation amounts to at most 
about ±10%. Moreover, the scale stability at small x does not significantly 
depend on tha. steepness of the gluon distribution. Consequently, the NLO 
results of ref. u seem to provide rather sound a theoretical foundation for a 
small-x gluon determination at HERA, despite the large total cm. energy 
s 3> 4m 2 which might suggest a destabilizing importance of ln[s/(4m^)] terms. 
At large x, x ss 0.1, on the other hand, the scale dependence of F 2 C is rather 
strong, especially at low Q 2 , presumably due to the large small- ij threshold 
contributions mentioned above. However, F£ is small in this region. 

The next issue we investigate is the locality in y of the gluon determination 
via F£ . The contribution from initial-parton momenta smaller than y max to 
F2{x,Q 2 ), denoted Jpy F£(x, Q 2 , y ma x), is presented in Fig. 3 for the GRV 
parton distributionsa. At scales /i rs JQ 2 + 4mj?, about 80% of F 2 C originates 
in the region y m i n = aa; ^ 3y m i n . Again the situation is very similar for the 
CTEQ2 parton densities Ej. Thus F 2 C allows for rather local a determination of 
g(y,[i 2 ~ \J Q 2 + Am 2 ). The partonic cm. energies in the region contributing 
80% to the structure function F£ are given by 77 < 5 (20), corresponding to 
s<60(180) GeV 2 , at Q 2 = 10(100) GeV 2 , respectively. This implies that for 
the Q 2 values under consideration here, the plateau region of and c' 1 ' at 
large 77 (cf. Fig. 1) does not play an important role. Resummation of large-?? 
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FE logarithms is thus not necessary, and not appropriate if it implies additional 
approximations in the more important small- s region fl. A similar observation 
has already been made inQ for /i ~ 2m c at low Q 2 . The latter scale choice leads 
however to a considerably wider important range of s at high Q 2 , somewhat 
in contrast to its motivation described above. 



1 


r F 2 ° (x 

• i- 


i 1 

Q 2 =10 QeV 


2 m 

; y 


<\ 

,x ) i 
i o -4 = x : 










10" 3 


10 " 






I 


io 2 : 


-2 

10 














jj. — 4m 2 i 
jj. = Q + 4m~ 






10 











= 1 

: f 2 c (x. 


Q 2 =100 GeV 2 ; y n,iiX ) 




1 . 


10"* = x 


l 






-i 

10 




Y 


10 " 


i — m 


2 = 4m 2 ' 1 




M 


2 = Q 2 + 4m 2 : 



LU LU LU 1U y " 1U 1U 1U 1U y 

Figure 3: The contribution of the initial-parton momentum region ax < y < j/ max to -F 2 C at 
small x for two choices of the scale /i, using the parton densities of ref.p. The arrows indicate 
the values of y max at which 80% of the complete results are reached for /i = •J Q 2 + 4m^ . 
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Figure 4: The x-deoendence of _F 2 C and F£ /F% at two fixed values of Q 2 , as expectedjrom the 
GRV gluon density!]. Also shown are F 2 C as obtained from the CTEQ 2MF parton setli2l and the 
charm mass dependence of the predictions, fi, = yj Q' 2 + 4m 2 has been employed. 

The expected absolute and relative sizes of F 2 C are displayed in Fig. 4. In 
contrast to the bottom contribution F 2 b which reaches at most 2 . . . 3%, F 2 C is 
large in the HERA smalLp region, making up up to a quarter of F2 as measured 
at HERA. See also ref. U. This size of F 2 C renders a reliable (fully massive) 
treatment mandatory in any precise analysis of Fi at small x. The sensitivity 
of F 2 C to the gluon density— is illustrated by the difference of the CTEQ 2MF 
(flat xg{x, [i 2 = 2.6 GeV 2 ) E3 and the GRV (steep gluon)l3 expectations. There 
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is quite some discriminative power of F 2 C , especially close to the lower end of 
the Q 2 range considered here, Q 2 « 10 GeV 2 . Moreover, by measuring up to 
about 100 GeV 2 in Q 2 , the rapid growth of yg(y <C 1,M 2 ) predicted by the 
Altarelli-Parisi equations can be rather directly tested down to y ~ 10~ 3 . A 
theoretical obstacle to an easy accurate gluon determination via the charm 
structure function is the dependence of F£ on the unknown precise value of 
the charm quark mass m c . A ±10% variation of m c also considered in Fig. 4 
leads to a ±15 . . . 25% (5 . . . 10%) effect at Q 2 = 10 (100) GeV 2 , respectively. 

To summarize: The NLO perturbative QCD approach to the charm struc- 
ture function F£ seems to be in good shape at small x where F£ is expected 
to be large (up to about a quarter of the total Fi in the HERA regime) , with 
scale variations of less than about ±10%. The situation is worse at large x, 
where the structure function is however small. F 2 C represents a clean gluonic 
observable in NLO and is well suited for rather local a gluon measurement at 
small- a;, with bulk of the result originating from gluon momenta within a fac- 
tor of three above the threshold value. Flavour excitation contributions from 
s ^ 4m 2 become important only at scales Q 2 higher than those relevant for 
small- x observations at HERA considered here. Despite the significant charm 
mass dependence of the results, especially at low Q 2 , useful constraints will be 
put on the proton's gluon density and its evolution at small x, if F 2 C data with 
an accuracy on the 10% level can be obtained at HERA. 
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